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I. INTRODUCTION 

With this paper we complete the K-matrix analysis 
of GAMS data on the reactions 7r~p ir^Tr^n M, rj-qn 
Hi rjrj'n [| which was started by the papers Ijfl]. The 
K-matrix analysis gives a rich information about meson 
states, thus helping the qq state classification and the 
search for exotic mesons. However for the restoration of 
the K-matrix amplitude, one needs to study a full set of 
open channels with sufficiently high statistics. It is the 
reason to include in our fit the data on 7r~p — > KKn |^] 
and ppiat rest) tt'^tt^'tt^, Tr^rji], Tr'^Tr^Ty 

K-matrix poles which are a subject of the present con- 
sideration differ from the amplitude poles in two points: 

(i) The states corresponding to the K-matrix poles do 
not contain components with real mesons which are in- 
herent in resonances. The absence of a cloud of real 
mesons allows one to refer conventionally to these states 
as the bare ones 

(ii) Due to the transition bare state — > real mesons — > 
bare state, the observed resonances are mixtures of the 
bare states. So, for the quark systematics, the bare states 
are primary objects rather than resonances. 

Coupling constants bare state — > real mesons are re- 
sponsible not only for the mixing of states but for the res- 
onance decays as well; the relations between couplings al- 
low one to restore the quark/gluon content of bare states 

The paper is organized as follows. 

In Sec. II we introduce a set of formulae which 
are used in the data fit. We present the S- and D- 
wave K-matrix amplitudes for the mass-on-shcU reactions 
TTTT — > TTTT, KK, TjT], rjrj' together with those for the mass- 
off-shell pion in the initial state: 7r7r(t) with t ^ m^. The 
K-matrix formulae for the final state interaction in the 
three-meson production process, pp three mesons, are 
presented as well. 

In Sec. Ill we write down the couplings for the tran- 
sition bare state — > two pseudoscalars, with the im- 
posed quark-combinatorics constraints both for gg-states 
(isoscalar and isovector) and for the glueball. Restora- 
tion of couplings in the fit allows us to determine the 
quark content of isoscalar states and to find a candidate 
for the glueball. 

Mesons which belong to the same qq nonet have ap- 
proximately equal masses; they also have approximately 
equal decay couplings. Besides, flavour wave functions 
for isoscalars of the same nonet are orthogonal. In Sec. 
IV we present the results of the fit with the imposed 
nonet-classification constraints. Fit of the 00+"'" wave 
confirms the result of Ref. while for the 02+"^", lO"'""*", 
12+"'" waves the K-matrix representation of amplitudes 
in the mass region below 1900 MeV is done for the first 
time. The restored bare states, together with those found 
in the Kn S-wave K-matrix analysis |pT|] , allow us to con- 
struct the l^Poqq and \^P2qq nonets unambigously, and 
for the 2^Poqq nonet two variants are possible which dif- 



fer in the mass of the lightest scalar/isoscalar state. 

The origin of the hghtest scalars, /o(980) and ao(980), 
is crucial for the nonet classification. These states are 
located near the KK threshold and give rise to the ques- 
tion whether these states are hadronic KK molecules. In 
Sec. 5 we present arguments based on the direct GAMS 
measurements together with the results of the performed 
K-matrix fit that the bare states from which /o (980) and 
ao(980) originate have qq nature. 

Short summary is given in Sec. 6. 

II. EXPERIMENTAL DATA AND K-MATRIX 
AMPLITUDE 

Here we briefly introduce the fitted data and set out 
the K-matrix formulae used for the data analysis. 

A. Experimental data 

Simultaneous analysis of meson spectra in the chan- 
nels IJP'^ = 00+"^, 10++, 02++ and 12++ is performed 
on the basis of the following data set: 

(1) GAMS data on the S-wave two- meson production in 
the reactions np — > n'^Tr'^n, rjrjn and rjrj'n at small nu- 
cleon momenta transferred, \t\ < 0.2 (GeV/c)^ IHI]; 

(2) GAMS data on the tttt S-wave production in the 
reaction np — > Tr'^n'^n at large momenta transferred, 
0.30 < \t\ < 1.0 (GeV/c)2 Q; 

(3) GAMS data on the tttt D-wave production in the re- 
action irp Tr'^TT^n, at small and large < \t\ < 0.5 
(GeV/c)2 I; 

(4) BNL data on TriV -> KKN §; 

(5) Crystal Barrel data on pp {at rest) tt'^tt^tt^, tt'^tt^t;, 
TT°riri 111; 



B. K-matrix amplitude and analyticity 

The K-matrix technique is used for the description of 
the two-meson coupled channels: 

A = K{I -ipK)-\ (1) 

where K is n x n matrix, where n is the number of 
channels under consideration and / is unit matrix. The 
phase space matrix is diagonal: pab — ^abPa- The phase 
space factor pa is responsible for the threshold singulari- 
ties of the amplitude: to keep the amplitude analytic in 
the physical region under consideration we use analytic 
continuation for pa below threshold. For example, the 
rjrj phase space factor pa ~ {1 — Am^/s)'^/'^ is equal to 
i{Am'^/s — 1)^/^ below the 7777 threshold (s is two-meson 
invariant energy squared). The phase space factors we 
use lead to false kinematic singularities at s = (in all 
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factors) and at s = (m^/ — m^)^ (in the 7777' space factor), 
but these false singularities which are standard for the K- 
matrix approach are rather distant from the investigated 
physical region. 

For the multimeson phase volume in the isoscalar sec- 
tor, we use the four-pion phase space defined as either 
pp or a (J phase space, where a denotes the S-wave tttt 
amplitude below 1.2 GeV. The result does not practicaly 
depend on whether we use pp or aa state for the descrip- 
tion of multimeson channel: below we provide the for- 
mulae and the values of the obtained parameters for the 
pp case, for which the fitted expressions are less cumber- 
some. The multimeson phase space in the sector / = 1 is 
taken in the form which, in its low-energy part, simulates 
the flop phase space. 



C. Isoscalar /scalar, 00^^, partial wave 

For the S-wave interaction in the isoscalar sector, we 
use the parametrization similar to that of Ref. : 

00 1 \ 



So 



s - mil 2 



So 



(2) 



where K^^ is a 5x5 matrix (a, 6 — 1,2,3,4,5), with the 
following notations for meson states: 1 — tttt, 2 = KK, 
3 = 7777, 4 = 7777' and 5 = multimeson states (four-pion 

state mainly at ^/s < 1.6 GeV). The gi""* is a couphng 
constant of the bare state a to the meson channel; the 
parameters fab and sq describe the smooth part of the 
K-matrix elements (sq > 1.5 GeV^). We use the factor 
(s — m1/2)/s to suppress the influence of the false kine- 
matic singularity at s = in the amplitude near the tttt 
threshold. 

The phase space matrix elements are equal to: 



Pais) 



a = 1,2,3, 



(3) 



where mi — m-^, 7712 — m^, 777,3 = "^r;, and 

p (s) = I at s > {mr, - mrj')^ 
' \ pi2 at s < (m^i — 771^')^ ' 



P41 = 
P42 = 0. 

The multimeson phase space factor is defined as 



(4) 



P5(s) 



P51 at s < 1 GcV 
P52 at s > 1 GeV^ 



P51 = Pa 



dsi f ds 



X M2r(si)r(s2) V(^+^7^^^F^4^ 

X [iM'-s2r+M'r\s2)r\ 

P52 = 1 . (5) 

Here si and S2 are the two-pion energies squared, M is 
/9-meson mass and r(s) is its energy-dependent width, 
r(s) — 7p5(s). The factor pQ provides the continuity of 
P5{s) at s = 1 GeV^. 

The following formulae describe tttt, 7777 and 7777' pro- 
duction amplitudes due to pion t-channel exchange: 

X (1 - ipA')^^^ , 6 = TTTT, ?7?7, 7777' , 

, , IGcV^ + so^ 



r(t),a 



fait) 



So 



X s - 



(6) 



ml/2) A 

Here TV is a normalization factor, Fpf{t) is the nucleon 
form factor, and D{t) is the pion propagator: 
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where A's, g' and /' are the fitted parameters. 



D. Isoscalar /tensor, 02^"'", partial wave 

The D-wave interaction in the isoscalar sector is 
parametrized by the 4x4 K-matrix where 1 — tttt, 2 — 
KK, Z ~ rjrj and 4 ~ multimeson states: 

K'^is)-Dais)^^''' 



Ml-s 

\ a. " 

lGeV2+£2\ n ( ^ 

+ /a6 \ Fibis) 



(8) 



Factor Dais) stands for the D-wave centrifugal barrier. 
We take this factor in the following form: 

^.(s) = ^^, a =1,2, (9) 



where ka = \J sj^ — ma is the momentum of the decay- 
ing meson in the cm. frame of the resonance. For the 
multimeson decay, the factor 1)5(5) is taken to be equal 
to 1. The phase space factors used are the same as those 
for the isoscalar S-wave channel. 
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E. Isovector /scalar, 10^"'', and isovector /tensor, 
12+^, partial waves 

For the amplitude in the isovector/scalar and isovec- 
tor/tensor channels, we use the 4x4 K-matrix with 1 
irrj, 2 — KK, 3 — irrj' and 4 — multimeson states: 



(a) (q) 

9a 'gl 



Kliis)^D.is)\y: 

\ a. " 

1.5 GeV^ + si 



+ fab 



Si 



Dbis) . 



(10) 



Here J — 0,2; the factors Da{s) are equal to 1 for the 
10"*""*" amplitude, while for the D-wave partial amplitude 
the factor Da{s) is taken in the form: 



Dais) 



kl + 3/r 



2 , a = 1, 2, 3, 



Diis) = 1 



(11) 



The elements of the phase space matrix in the isovector 
sector are defined as 



Piis) = 



Pii 

Pl2 



at 
at 



s > (rrirf — m-^Y 
s < (rrirf — rriTrY ' 



pi 

Pl2 = 0, 



P2is) 



(12) 



(13) 



Psis) = 



P31 
P32 



at s > (mjj — mTr)^ 
at s < (m'^ — m^)^ ' 




The multimeson phase space factor pi{s) is taken in the 
form which simulates the pao phase space factor below 
s = 2.25 GeV^: 



Pi{s) 



P41 at {nirj + 3m^)2 < s < 2.25 GeV^ 



P42 
P43 



at s > 2.25 GeV 

at s < {rriri + 3m^)^ 



PAl = 



1 — {rriri + im-^Y / s 
1 - (m^ + 3m^)2/2.25 GeV^ 



5/2 



P42 = 1 , 
P43 = . 



(15) 



F. Three-meson production amplitudes 

The amplitudes pp {at rest) — s- tt'^tt'^tt'^, Tr'^T/Ty which 
correspond to the production of the two-meson isoscalar 
states are equal to: 



ICC mesons 



^i(23) + A2(13) + A3(12) , 



where the amplitude Ak{sij) stands for diagrams with an 
interaction of particles in the intermediate states and the 
last interaction being of the particles i and j, while the 
particle A: is a spectator. We suppose, as in the previous 
papers [^J^ , that pp annihilates at rest from the ^S'o-level. 
The following form is used for the two-particle interaction 
block: 



Ai(23)= ^./(23) i^p%,„(s23) 



J=0,2 



^{^-WK^-'{^23))1 



-1 
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(16) 



Here h = 7r°7r° stands for 7r°7r°7r° production, and b = rjrj 
for ir'^rjr]. The centrifugal barrier factor Xj is equal to 
1 for the production of the S-wave resonance. For the 
D-wave resonance production, this factor is: 



^2(23) = i(3cos2 912-1) 



3/i?2' 



(17) 



where 9 12 is the angle between particles 1 and 2 in the 
rest frame of the particles 2 and 3, pi is the momentum 
of the particle 1 in this frame and R characterizes the 
annihilation radius. The X-matrices which describe the 
prompt resonance production in the pp annihilation have 
the following form: 



A^;-l[00]gi") 



^ Ml - s 



1 GeV^ 



So 



So 



-m^/2 



(18) 



is-]- 



f^pp 



p7T,a [0; 2j 



IGeV-'-hsoX 

X 7 -)Da{s^J). 19 

Sij + Sq J 

The TTTTTT production amplitude is completely described 
by eqs. (|l^)-(|lS|) because of the amplitude symmetry 
under the rotation of pion indices k. 

The part of the amplitude {at rest) —>■ mr'^Tr'^, which 
corresponds to the production of isoscalar resonances, 
reads: 

Ai(23)= ^ Xj{23) k^J^Js23) 



J=0,2 



{1-ZPK"'{S23))1 



. -1 

2fc 



(20) 
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where 



£>oo / \ ^ Appi [OOlsi '' 



Ml - 



1 GeV^ + So 

Si-j + So 



(21) 



Parameters App^[OJ] and (f)pp-„\fiJ] (or App^[OJ] and 
0ppi)[OJ]) may be complex magnitudes with different 
phases due to the three particle interactions. 

The part of the amplitude, which corresponds to the 
production of the isovector resonances in the reaction 
pp [at rest) —>■ rjrjTr'^ , is written as Ai{23) + ^2(13) and 



^2(13) = ^•'(13) K'J,,aisi3) 



J=0,2 



where 



rjTT 



(22) 



~\~(f>ppri,a \^J\ 



1 GeV-^ + si 

Si] + Sl 



Dais.j). (23) 



The production of isovector resonances in the reaction 
pp {at rest) -K^-K^rj has the form Ai(23) + ^2(13) and 



^2(13)= ^J(13) i^pY..a(si3) 



,7=0,1,2 



X (l-ipi^^^(si3)^, , 5 = r;7r", (24) 
where K^p^ ^ is given by Eq. (|2^) with the replacement 

^p^ill'^] ^ App^ril'^] and 4>ppn,a[^J] ^ 0pp7r,a[lJ]- 




FIG. 1. Quark-antiquark loop diagram which determines 
the glueball width (a); diagrams for the decay of a gg- meson 
(b) and a glueball (c),(d) into two gg- meson states. 

We calculate the ratios of the decay coupling constants 
in the framework of the quark combinatoric rules which 
were previously suggested for the high energy hadron pro- 
duction and then extended for hadronic J/^' decays 
]l5[ . Calculations of the decay coupling constants for the 
glueball and isoscalar/scalar gg-mesons were performed 
in Refs. [p|,p|jlO|]. The decay couplings for isoscalar and 
isovector mesons are given in Tables | and ||. 

Isoscalar meson decay couplings depend on the non- 
strange/strange component ratio of the decaying meson 
given by the mixing angle <&: 



V', 



flavour 

decaying meson 



nnCos^ + ssSin^ 



(26) 



Here nh = {uu + dd) / ^/2. It allows us to restore $ and 
at the same time to determine the decay couplings. 

The glueball decay couplings obey the same ratios as 
the isoscalar/scalar gg-meson couplings with the mixing 
angle equal to: 



III. QUARK-COMBINATORIC RULES FOR THE 
DECAY COUPLINGS AND THE qq CONTENT 
OF MESONS 

The decay couplings of the gg-meson and glueball to 
the two mesons are determined by the diagrams with qq- 
pairs produced by gluons. Figs. |^(b),(c) provide an exam- 
ple of diagrams which contribute to the leading terms in 
the 1 /N expansion and Fig. |l|(d) is an example of dia- 
grams for the next-to-leading contribution. The produc- 
tion of soft qq pairs by gluons violates flavour symmetry, 
with the following ratios of the production probabilities: 



$ = $ 



glueball 't^'^^gluebaU 



(27) 



uu : dc? : ss = 1 : 1 : A 



and A = 0.4 - 0.8 131. In our fit we fix A = 0.6. 



(25) 



It follows from the two-stage decay of the glueball , 
see Fig. |^(c): an intermediate qq-sidXe in the glueball 
decay is a mixture of nn and ss quarks produced in pro- 
portion given by Eq. (H). We fix ^giuebaii = 25° ± 5°. 

The coincidence of the glueball decay couplings with 
those of qq meson at $ = ^giuebaii points out that there is 
no simple signature of a determination of glueball state: 
searching for the glueball, we should perform a full qq- 
classification of mesons, thus an existence of an extra 
state for the gg-classification is an indication of the ex- 
otics. 

The normalization in Table | is done in such a way 
that for the glueball decay the summation of couplings 
squared over all channels is proportional to the probabil- 
ity to produce a two-quark pair, (2 + A)^, see Eq. (p^). 
So, 
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channels 

E 5^(c)/(c) = ig^(2 + A)^ 



(28) 



channels 



Here /(c) is the identity factor and c = tt^tt^, tt+tt", 
ii'+ii'^, and so on (see Table |). With this normahza- 
tion gc/G ~ l/N^. The experience of the quark-gluon 
diagram calculations teaches us that the factor 1 /Nc ac- 
tually leads to a suppression of the order of 1/10: in the 
fitting procedure we impose a restriction \gG/G\ < 1/3. 

The nonet classification of isoscalar mesons is based on 
the following two constraints: 

(1) The angle difference between isoscalar nonet part- 
ners should be 90°. For this value the corridor ±5° 
is allowed in our analysis: 



$(1) - $(2) = 90° ±5° . 



(29) 



(2) Coupling constants g of Tables | and |l| should be ap- 
proximately equal to each other for all nonet part- 
ners: 



g[fAl)] c g[fA2)] c g[aj] c g[Kj] 



(30) 



The conventional quark model requires exact coincidence 
of the couplings g but the energy dependence of the loop 
diagram of Fig. |l](a), B(s), may violate this couphng 
constant equality because of the mass splitting inside a 
nonet. The K-matrix coupling constant contains an addi- 
tional s-dependent factor as compared to the coupling of 
the N/D-amplitude 0: g^{K) = g^{N/ D)/{1 + B' {s)). 
The factor (1 + B'{s))~^ mostly affects the low-s re- 
gion due to the threshold and left-hand side singularities 
of the partial amplitude. Therefore, the coupling con- 
stant equality is mostly violated for the lightest 00++ 
nonet, I^Pq QQ- We allow for the members of this nonet 
1 < 5[/o(l)]/5[/o(2)] < 1.5. For the 2^Po qq nonet mem- 
bers, we put the two-meson couplings equal both for 
isoscalar and isovector mesons. The equality of coupling 
constants is also imposed for tensor resonances. 



IV. DESCRIPTION OF DATA AND THE 
RESULTS 

The performed K-matrix fit gives a good description 
of the data, see Figs. |^|[ The values for the fit are 
shown in Table [II and parameters of the fit are presented 
in Tables IV-VII. Below we single out the main results 
of our fit. 



A. IJ^^ = 00++ wave 

The present fit confirms the results obtained in the pre- 
vious analysis of the 00++ wave . For the description 



of the 00++ wave in the mass region below 2000 MeV, 
five K-matrix poles are needed (a four-pole amplitude 
fails to describe well the data set under consideration). 
Correspondingly, five bare states are found: 

^bare^720± 100), 
^fiavour ^ 45 _^ 0.1)wl - (0.89 ± 0.05)ss , 

/obare(i230±50), 
^fiavour ^ (^o,g+of)nn + (0.45t^:?)ss , 

/obare(i260±30), 
^fiavour ^ (0.93+0-02)^^ + (0.37tg:gg)ss , 

/□^^"^^(leooiso), 

^fiavour ^ (0 95 ^ q q^)^- ^ (0.3^014)^- ^ 

/obare(i8io±30), 

(0.10 ± 0.05)ml -I- (0.995lg;^?^)ss, 

Solution /), 
(0.67 ± 0.08)nn - (0.74 ± 0.08)ss, 
(Solution //). 



fiavour 



(31) 



Experimental data used in the fit do not fix unambigously 
the fiavour wave function of /|^^'^*^(1810 ± 30): two solu- 
tions are found for it. 

The scattering amplitude has five poles in the energy 
complex plane, four of them correspond to relatively nar- 
row resonances while the fifth resonance is very broad: 



/o(980) 
/o(1300) 
/o(1500) 
/o(1530) 

/o(1780) 



(1015 ± 15) - i(43±8) 
(1300±20)-i(120±20) 

(1499 ±8) - i(65± 10) 

(1530j;^^o) - i(560 ± 140) 

(1780 ± 30) - i{UO ± 20) MeV, 

(Solution /), 
(1780 ± 50) - i(220 ± 50) MeV, 
(Solution //). 



MeV, 
MeV, 
MeV, 
MeV, 

(32) 



The broad resonance is crucial for the description of the 
00++ wave, being responsible for large interference ef- 
fects which are seen in different reactions. Namely, the 
resonance /o(980) reveals itself as a dip in the S-wave 
TTTT TTTT-spectum, Fig. H(a), and as a sharp peak in the 
TTiT{t) nn spectra at large \t\, Fig. |[ The resonance 
/o(1300) is seen in the 7r7r(i) — > tttt spectra at large \t\ 
as a well shaped bump. Fig. ^, while in the tttt tttt 
and TTTT — > KK spectra it reveals itself as a shoulder, 
Figs. H and ||. /o(1500) is seen as a dip in the tttt tttt 
and TTTT rji] spectra. Figs. ^, ||, and as a peak in 
pp{at rest) Tr^Tr^Tr" reaction. Fig. ^. In all these ap- 
pearances of /o(980), /o(1300) and /o(1500), their inter- 
ference with /o(1530+25o) plays a decisive role. In the 
case of large interference effects it is useful to display the 
amplitude on the Argand-plot. The Argand-plots for the 
amplitudes tttt —^ tttt, tttt — > 7777, tttt — > KK, tttt — > rjr]' 
and inr{t) — > tttt are shown in Figs. ^ and ||. 
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Four bare states of Eq. (^ij) can be naturally classified 
as nonet partners of the qq multiplets I^Pq and 
The fifth bare state is superfluous for the qq classification 
being a good candidate for the lightest scalar glueball. 
Eq. ( |3l] ) gives two variants for the glueball: either it is a 
bare state with mass near 1250 MeV or it is located near 
1600 MeV. Correspondingly, after having imposed the 
constrains ( p9|) and (30), we found the following variants 
of the nonet classification. For the solution I: 

I. fl>^^^{720) and f^^^^{1260) are I^Pq nonet part- 
ners, 

/^^'^'^(leOO) and fl'^'^^ilSlO) are 2^Po nonet part- 
ners. 



/obare(i230) is a gluebaU. 



Within solution II, two variants describe well the data 



set: 



II-l. /,J'^i"^(720) and fl'^^'^ (1260) are l^Po nonet part- 
ners, 

jbare(;LgOO) and f^'^'"'^ (1810) are 2^Po nonet part- 
ners, 

/obare(i230) is a gluebaU; 



II-2. /obare(720) and f)^^^^{1260) are l^Pa nonet part- 



ners. 



jbare(j^230) and f^^^^{1810) are 23^0 nonet part- 



ners. 



f^^^^{1600) is a gluebaU. 



Tables [V-^ present parameters which correspond to 
these three variants. 

Lattice calculations of the gluodynamic glueball 
give the mass of the lightest scalar state in the re- 
gion 1550-1750 MeV that coincides with the variant II- 
2. However, it should be emphasized that the state 
jbare(;LgOo) can not be identified as a pure gluodynamic 

glueball because fo'^^^'s contain the gq-components re- 
lated to real parts of the loop transition diagrams: this 
problem is discussed in detail in Rcfs. |l^jl^,^. An ex- 
traction of the gg-component from f}j'^^'^ {1600) leads to 
the mass shift of the state which is not large according 
toRefs. MM: /obare(i600) ^ /o^""" (1633). 



1520±40 MeV. A similar result is obtained in the present 
K-matrix approach. 

For the description of the isovector/isoscalar scattering 
amplitude, we use the two-pole 4x4 K-matrix with two- 
meson coupling constants given in Table ||. 

In the first stage of the fit, the coupling of the lightest 
ao-state was allowed to vary in the interval bounded by 
g[jbare^720)] and .g[/o^arc(;^260)]. In all the variants of 
the fit, the two- meson coupling constant of the lightest 
state, g[aQ^^^ (lightest state)], appeared to be very close 

to the couphng constant g[fl'^^'^ {1260)]: in the final fit, 
in line with the constraint of Eq. (30), we fix these cou- 
plings equal to each other. The two-meson coupling of 
the next isovector/scalar is fixed to be equal to the cou- 
plings of the 2'^Pq isoscalar/scalar states. 

The fit gives two solutions for the 10++ wave 
which practically coincide in terms related to the 
resonance/bare-state sector and differ in background 
terms. Parameters for the 10++ wave and the pole po- 
sition are given in Table VI. We have for the resonance 



positions and the bare states, correspondingly: 

ao(980) 
ao(1450) 



±6)-i(46±10) MeV, 
(1535±30) -i(146±20) MeV (33) 



and 



a^'^^^(964± 16) 



,bare 



(1670 ±80) 



(34) 



But these two solutions give different predictions for the 
scattering amplitudes: for the first solution (without K- 
matrix background terms) the ttt] — > m] scattering am- 
plitude squared (see Fig. ^ has a dip in the region 
of ao(1450) due to the destructive resonance interference 
with the background, while for the second solution (with 
the K-matrix background terms) a dip appears at 1100 
MeV. In the present fit, the information on the isovec- 
tor/scalar wave comes from Crystal Barrel data only. 
These data being rather sensitive to the pole structure 
provide poor information about K-matrix background 
terms: it is a source of ambiguities in our K-matrix so- 
lution. But, let us stess, the description of other partial 
waves practically does not depend on whether the first 
or second solution is used: the variation of parameters is 
within the given errors. 



B. IJ^^ = 10++ wave 

Two isovector / scalar resonances are well seen in the pp 
annihilation into three mesons p|,|6|,p|, pj|j2C| ] . The lightest 
one is the well known ao(980), while the next resonance 
is the newly discovered ao(1450) with mass 1450 ± 40 
MeV and width F = 270 ± 40 MeV, as is given by the 
Particle Data Group ||2^. Let us note that in fitting 
the last high statistic Crystal Barrel data with the T- 
matrix method used for this wave |^,^,^ the mass of 
this resonance appeared to be a bit higher and equal to 



/J^C' ^ 12++ 



Similar to the isovector/scalar case, the 4x4 two-pole 
K-matrix is used for the description of the 12++-wave. 
Coupling constants of bare states and t he p oles of the 
scattering amplitude are given in Table VIl. We have 
determined two bare states: 



^bare/ 



'(1314 ±7), a|'^^^ (1670 ±75). 
The poles of the amplitude are located at 



(35) 
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a2(1320) 
02(1640) 



> (1309 ±6) 
(1640 ±50) 



i(58 ± 6) MeV, 
i(122± 18) MeV 



(36) 



The lightest state is a well known 02(1320) resonance, 
with mass 1318 ± 1 MeV and width T = 107 ± 5 MeV, 
according to Ref. plj . 

When fitting Crystal Barrel data on the reaction 
pp{at rest) — > titttt, the introduction of the isovec- 
tor/tensor resonance with mass 1600-1700 MeV makes 
an appreciable improvement of the Dalitz plot descrip- 
tion in this region. 



D. IJ^^ = 02++ wave 

The two lightest isoscalar/tensor states, /2(1270) and 
/2(1525), are well known: they are members of the nonet 
l^P2qq. Crystal Barrel data point out the existence 
of the resonance /2(1565) which helps to describe the 
pp — > TT^TT^TT^ Dalitz plot lu the region of large two-pion 
masses [p|, p7| , p0[ . Because of that, we also begin our anal- 
ysis introducing a three-pole K-matrix amplitude. How- 
ever, after imposing the nonet constraints on the 1^P2 
states, see eqs. (28) and (29), we obtain that the cou- 
plings of the third state turn out to be negligibly small. 
Although the description of the reaction pp — > ij-qir'^ be- 
comes a bit worse under the imposed constraints (about 
0.1 per degree of freedom for x^), the description of the 
react ion pp 7r°7r°7r° (where /2(1560) is seen according 



to p, |17y2(]| ]) improves x by 0.07 giving practically the 
same total x^- 

/2(1560) is not seen in GAMS data; that gives a strong 
restriction on the partial width of the resonance decay 
into TTTT channel, less than 20 MeV. 

In our final fit, we have used the two-pole K-matrix 
amplitude with the nonet constraints; parameters for this 
fit are presented in Table VII . 

The K-matrix fit gives the following bare isoscalar / ten- 
sor states, the members of the ■^P2 nonet: 



/|^^'^e(1235± 10), /|^^i"6(1530± 10), 
$[/|^^'^'^(1530)] = 86° ±5° . 



(37) 



The K-matrix 02+^ 
mass values: 



amplitude has poles at the complex 



/2(1270) -> (1262 ±6) 
/2(1525) ^ (1518 ±9) 



-i(90±7)MeV, 
i{71 ± 10)MeV. 



(38) 



These values should be compared with masses and half- 
widths of Particle Data Group |2^ which are, corre- 
spondingly: (1275 ± 5) MeV, (92.5 ± 10) MeV and 
(1525 ± 5) MeV, (38 ± 5) MeV. The width of /2(1525) 
found in our fit appears to be much larger than one 
given in PDG. It is quite possible that in fitting the 
present data set we cannot resolve a possible D-wave dou- 
ble pole structure in the region of 1530 MeV caused by 



the /2(1525) and /2(1560) resonances, for they are lo- 
cated near the edge of the phase space for Crystal Barrel 
data, while GAMS data give a restriction only on the cou- 
plings to TTTT channel. It is possible that the additional 
information from Crystal Barrel data on KKn produc- 
tion together with GAMS H] and VES data 
production will clarify this point. 



E. Nonet classification 

The results of the performed analysis together with 
the results of the K-matrix analysis of the Ktt S-wave 
|]lT| allow us to construct the lightest scalar qq nonet 
uniquely as 

l'Vo:/obare(720±100), 
/obare(i260±30), 
abare(960±30), 

i^bare(i220t^°o), 
$[/obare(720)] 



-70° +^° 



(39) 



The hghtest scalar, fl'^^'^{720 ± 100), is dominantly a 
ss state with mixing angle close to the ideal octet one, 
<I>idoai octet = —55°. The situation with the lightest scalar 
nonet is similar to that with the lightest pseudoscalar 
nonet, where the mixing angle for the 77-meson is also 
close to the 'I'idoai octet: this definitely indicates the de- 
generation of the hghtest 00++ and 00~+ states. 

The multiplet of the lightest tensor states appears as 



l'P2 



ybare(i240±10), 
/2bare(i522±io), 
abare(i3ii±3)^ 

ii:2*(1430) 
$[/2^a'^'5(1240)] = -10° ±3° 



(40) 



The K-matrix analysis of the irK D-wave is not done: 
the (J = 2)TrK resonance with mass 1431 ± 3 is reported 
in Ref. ; we have used this resonance to complete the 
multiplet (40). 

Our analysis gives two variants for the 2 ^Pq gg-nonet: 
First variant: 



2^Po : /(J'^^'*^(1600±50), 
/obare(i8io±30), 
a^^^^(1650±50). 



ii:obare(i885+50 



lOOJ 



$[/obare(^810)] = 84° ± 5° 



(41) 



The state k}^^"^^ {ISSbtl^o) is fixed by the analysis 0] 
of the Ktt 5-wave. In this variant the lightest glueball 
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state is fl'^'^^ {1230^11°)- In the second variant the light- 
est glueball state is identified as /(^*'-^(1600), namely: 



/obare(i8io±30), 

aj5'are(i650±50), 
Kbare(i885+50^)^ 

$[/obare(i8io)] = 44°± 10° 



(42) 



V. RESONANCES /o(980) AND ao(980) : ARE 
THEY KK MOLECULES? 

First, let us discuss the origin of /o(980). GAMS data 
for the /o(980) production at large \t\, see Fig.|3| directly 
demonstrate that this resonance has a hard component, 
while the location of the pole near the KK threshold 
definitely says that its kaon component is a long-range 
one. The existence of the long-range component gives 
rise to discussion about molecular structure for this state 
p5| . The problem to discuss is how substantial are these 
components in the formation of the resonance. Remem- 
ber that the short-range component (with r < 1 fm) is 
a subject of quark/gluon considerations and quark sys- 
tematics. 

The resonance /o(980) corresponds to the two poles 
located at (in MeV units): 

AI = 1015 — i46 (// sheet, under tttt — cut) , 
M = 936 - i238(/// sheet, under tttt and KK cuts). 

The second pole appears due to well-known double-pole 
structure caused by the KK -threshold (see, for exam- 
ple, 1^]), while the first pole at M = 1015 - i46 MeV 
generates the leading irregularities in tttt spectra. 

The restored K-matrix amplitude allows one to see the 
role of the KK component in the formation of /o(980), 
thus clarifying if this resonance is a descendant of a qq 
state or is a molecular-type system. To this aim, let us 
switch off the /o(980) decay processes (transitions into 
TTTT and KK) and look at the dynamics of pole posi- 
tions, with gradual onset of couplings. For the gradual 
change of couplings we performed the replacement in the 
K-matrix OO+^-amplitude: 



(a) 



Parameter ^ varies in the interval: 
< ^ < 1 . 



(43) 



(44) 



At ^ ^ the decay channels are switched off, and we 
have a bare state, while at ^ = 1 the real case is restored. 

At ^ ^ the masses of the lightest scalar bare states 
are 650 and 1260 MeV (the positions of the K-matrix 



poles). The trajectories of states with increasing f are 
shown in Fig. |l^. 

The crucial point is what component, tttt or KK, is 
mainly responsible for the mass shift from 650 MeV to 
1020 — i48 MeV. We can elucidate this point, switching 
off the KK component and leaving tttt untouched, and 
vice versa. In the first case the mass of /o(980) state is: 



M (without KK) = 974 - illb MeV . 

One sees that the mass shift 

^KK M(C = 1) - Af (without KK) = 
41 + i67 MeV 



(45) 



(46) 



is not large: the ii'^-component which is responsible for 
the value of S does not play an important role in the for- 
mation of the mass of /o(980). In the second case, when 
the TTTT component is switched off, we get the nearest state 
to the KK threshold, which is located at: 



M(without tttt) = 810 - ilO McV 
So the mass shift is 



(47) 



= 1) - A/(without tttt) = 

205 + m MeV, (48) 

being much larger than (5^^. The transition into real 
pions. 



/bare(720) -> ^ 



(49) 



is mainly in charge of the mixing of the lightest scalar qq 
state with other scalars thus shifting its mass by chance 
to the region of the next threshold, KK. The KK com- 
ponent of /o(980) is of the molecule- type: relative kaon 
momenta are small, so relative distances are large. But, 
let us stress again, the two-kaon component does not play 
a crucial role in the formation of the mass of /o(980). 



VI. CONCLUSION 

We have performed the K-matrix anlysis of GAMS 
data on the S- and D-wave tt^tt^, 7777 and 7777' data, to- 
gether with data obtained by BNL and Crystal Barrel 
Collaboration. Partial amplitudes for the states OO"*""*", 
02++, 10++ and 12++ are investigated in the mass region 
up to 2000 Me V, th e poles of these amplitudes are found, 
see Tables IV- VII. Pole terms of the K-matrix are re- 



stored, i.e. the bare states with quantum numbers 00++, 
02++, 10++ and 12++ are found. The quark content of 
these bare states is determined, based on the relations 
between coupling constants of the decays: this allows to 
restore the quark nonets I^Pq, and 1^P2- The per- 
formed analysis confirms the result of Ref. [6] which is 
based on the K-matrix analysis of the 00++ wave only: in 
the region 1200-1600 MeV there exists a scalar/isoscalar 
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state which is superfluous for the qq systematics. This 
state is a good candidate for the Hghtest scalar glueball. 

The analysis indicates the degeneration of the lightest 
00++ and 00"+ states. 
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FIG. 2. 

The TTTT — !■ TTTT S'-wRvc amplitude module squared |l| (a), the D-wave amplitude module squared (b), SD-correlator 
(c) and the phase difference between S and D- waves (d); the events are collected at the momentum transfer squared 
\t\ < 0.20 GeVVc' [|. 

Solid curve corresponds to solution II-2 and dashed one to solution I. 
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FIG. 3. 

Event numbers versus invariant mass of the TTTr-system in the S-wave for different t-intervals in the tt~p 
reaction Q. SoUd curves correspond to solution II-2 and the dashed curves to solution I. 
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FIG. 4. 

Event numbers versus invariant mass of the TTTr-system in the D-wave for different t-intervals in the tt~p —^ Tr'^ir'^n 
reaction The sohd curves correspond to solution II-2 and dashed one to solution I. 
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FIG. 5. 

The 5-wave amplitudes squared for transitions: (a) tttt KK (b) tttt ^ 7777 Q and (c) tttt rjrj' [^. The sohd 
curve corresponds to solution II-2 and the dashed curve to solution I. 
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FIG. 6. 

Mass projections of the Dalitz plot on the two-meson invariant mass for Crystal Barrel data. The curve corresponds 
to solution II-2. 
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FIG. 7. 

Argand plots for the S-wave scattering amplitudes in solution II-2: tttt rnr (a), tttt KK (b), tttt rjij (c) and 
WW — > 777/' (d). 
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Argand plots for the Tnr{t) tttt S-wave S-wave scattering amplitudes at different t. 
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FIG. 9. 

The squared S-wave nr] nr) scattering amplitude: solutions 1 (a) and 2 (b) for the nrj nr) scattering amplitude. 
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FIG. 10. 

Location of 00++ amplitude poles in the coniplcx-\/s plane {M = Re^/s, —T/2 — Im^) after replacing (/" ^9ai 
on the sheet under the tttt cut (a) and on the sheet under tttt and KK cuts (b). The case ^ ^ gives the positions 
of masses of bare states; ^ = 1 corresponds to the real case. The point 1 corresponds to ^ = 0.4, 2 to ^ = 0.6 and 3 
to ^ = 0.9. 
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TABLE I. 

Coupling constants given by quark combinatorics for a gg-meson decaying into two pseudoscalar mesons in the 
leading terms of the 1/N expansion and for glucball decay in the nc;xt-to-loading terms of the 1/N expansion. $ 
is the mixing angle for nn and ss states, and is the mixing angle for r] — rj' mesons: r} = nn cos Q — ss sin 6 
and rj' = nnsinO + sscosS. Glueball decay couplings in the leading terms of 1/A'' expansion are obtained by the 

rei)];u-omout.s /y/\/2 cos (I> — . ,f/siu<I> — \^G^ . 



Ch&imel 


The ^g-meson decay 

couplings in the 
1 pari 1 TIP" t.prms of 1 / AT 

expansion 


Glueball decay couplings 
in the next-to- 
leading terms of 1/iV 
expansion 


Identity 
factor in 


TT TT 


g cos$/-s/2 







1/2 




g cos$/%/2 







1 


K+K- 


^(V^ sin $ + cos *) /VS 







1 


K°K° 


5(1/2 sin $ + \/A cos *) /VS 







1 


nn 


g (cos^ e cos $/a/2 + a/a sin $ sin^ 0) 


2.c/G(cose - V 


1 sin 9)2 


1/2 


nn' 


g sin 6 cos 6 (cos 3>/\/2 — %/A sin 3>) 


2gG (cos e - sin @) (sin 9 + ^/^ cos Q) 


1 


n'n' 


g (sin^ e cos $/v^ + sin * cos^ 6) 


2gG{sm.Q + ^J 


A cos 6)2 


1/2 


TABLE IL 

Coupling constants given by quark combinatorics for scalar mesons and Og decaying into two pseudoscalar mesons 
in the leading terms of the 1 /N expansion. 


Channel 


The su- meson 
decay couplings 


Channel 


The du- meson 
decay couplings 






r?7r~ 




s^cose 




din 






5 ^sin9 


K-f) 


g ^(cosG - V2Asine) 


K°K- 




y 2 




g ■^(sine + A/2Acose) 
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TABLE III. 

^ values for the K-matrix solutions. 





Solution I 


Solution II-l 


solution II-2 


Number of points 


Crystal Barrel data Q 












1.52 


1.41 


1.42 


1338 


pp n'^r]!] 


1.57 


1.60 


1.59 


1798 


pp — > TT^TT^j; 


1.38 


1.43 


1.43 


1738 


+ - u u 

TV TT — > TT TT 










S-wave GAMS data 


1.47 


1.71 


1.59 


70 


D-wave GAMS data ^ 


1.63 


2.16 


2.14 


56 


SD-correlator jlj 


1.82 


2.26 


2.12 


47 


t-dependent GAMS data [§ 










0.00 < \t\ < 0.20 


3.03 


3.42 


3.37 


21 


0.30 < \t\ < 1.00 


2.64 


3.25 


2.98 


38 


0.35 < \t\ < 1.00 


1.30 


1.55 


1.44 


38 


0.40 < \t\ < 1.00 


2.75 


2.48 


2.79 


38 


0.45 < \t\ < 1.00 


1.92 


1.49 


1.67 


38 


0.50 < \t\ < 1.00 


2.29 


1.85 


2.04 


38 


GAMS data ^ 










7VTV ~* rirj 


0.70 


0.97 


0.87 


16 


TTTT rjri' 


0.49 


0.65 


0.64 


8 


nn — > KK 










BNL data § 


0.88 


0.77 


0.97 


35 



TABLE IV. 

Masses, coupling constants (in GeV) and mixing angles (in degrees) for the /,|^^'"^-resonances for solution I. The errors 
reflect the boundaries for a satisfactory description of the data. II sheet is under the tttt and 47r cuts; IV sheet is 
under the tttt, 47r, KK and rji] cuts; V sheet is under the tttt, An, KK, rjrj and rjrj' cuts. 





Solution I-l 




Q = 1 


a = 2 


Q = 3 


a = 4 


a = 5 


M 


0.65lt:J^« 


l-247l;-g 


1 2'S3+— 


1.684l;«° 


1. / 040 






0.5971:?^° 


0.8791°^° 


0.702l;°26o 


0.7021;°^° 


9G 





-0.1351:^0 











(a) 
95 





0.9441-0° 





0.8981 °™ 


0.3021;-° 




-(71.5tL) 


21. 5l^ 


U.ltf 


-6.011° 


891^5 




a — TTTT 


a = KK 


a = rirj 


a — rjrj' 


a = Itt 


/la 






0.5011:1°° 


0.448l;l°° 


-0 129+°''° 






-0 259+°-''^ 


b = 2,3,4,5 

= 


r, rjyr + O.lOO 

"-'•^'•^-0.100 


SO = 3.25l5"o 


II sheet 


imetZ 

-iiOMStZ) 




Pole position 






IV sheet 




1 QAQ+UiU 

-i{0.13StZ) 


1.496l-°« 
-^{0.059tZ) 


1.6701:1°° 
-z(0.760l:??°) 




V sheet 


1.7751-- 
-i(0.056l:°-,) 
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Masses, 
II-2. 


coupling constants (in 


GeV) and mixing an 


TABLE V. 

gles (in degrees) for the /g'^'-'^-resonances 


for solutions II- 1 and 




Solution II- 1 




a = 1 


a = 2 




Q = 3 


a = 4 


a = 5 


M 


0.65li;J^° 


l-246l;J^° 




1.2631:°!^ 


1.5951:°^° 


1.8321:°^° 


„(a) 




0.375l;°I° 




0.923l:H 


0.4241:°^,° 


0.4241:°^° 


9G 





-0.017t:°?° 













(a.) 

% 





0.705l;q 







0.552l:QyQ 


— 0.557l:o7o 




-(70.ll?5) 






18.3l^ 


20.6l?^ 


-64.411° 




a = TTTT 






a = ijrj 


a = rjrj' 


a = 47r 


/la 


0.4401:^ 


-0.0641:1°° 




0.3871:1°° 


0.4191:1°° 


-0.165l:°«° 






/6a=0 




6 = 2,3,4,5 










U.ZO»_Q 046 






-0.284l°:l°° 


So = 3.28lJ°o 


II sheet 


i.oi7t:r 

-z(0.049l:°°^) 






Pole position 






IV sneet 




i.3iii:-; 
-i(o.ii7i:;]^^) 




1.0UU_ 006 

-i(0.063l:°'«) 


1.41/ U_ 100 

-z(0.545l:-») 




V sheet 


1.8141;"- 
-i(0.082l;"o?^) 




Solution II-2 




a = 1 


a = 2 




a = 3 


a = 4 


a = 5 


M 


0.65ll:J^° 


1.219l:i^S 




1.2671:°!^ 


1.584l:°i° 


l-817l:°ofo 




1 oi:i+.100 

1.001_ J^QQ 


0.4351°^° 




o.9oii:°«° 


0.4331:155° 


0.4351:°™ 


go 













-0.0051:°^° 





(ol) 
96 





0.7191:1°° 







0.5421:°^° 


-0.5121:°™ 


$a 


-(69.5l?5) 


40.71^ 




19.61^° 


20.0l?^ 


-5411° 




a = TTTT 


a = KK 




a = rjri 


a = rjri' 


a = 47r 


/la 


0.459l;}°» 


0.0461:1°° 




0.4051:1°° 


0.4201:1°° 


-0 214+ °''° 






/6a =0 




6 = 2,3,4,5 










-0.24ll°:°« 






n O7q+0 l00 
"■^''J-O.lOO 


so = 3.05lJ°o 


II sheet 


i.o2ol:«°« 

-i(0.048l:r8) 






Pole position 






IV sheet 




1.304l:H 

-i(o.ii8l:°i^) 




l.dUd_ 008 

-i(0.063l:°°^) 


1.420l:i- 
-i(0.540l:°^°) 




V sheet 


i.809i-;;i';; 

-'"(0.080+:;;^;;) 
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TABLE VI. 

Masses and coupling constants (in GeV) for ao resonances. The star denotes that the parameter is fixed. 





ao-resonances 


without K-matrix background term 






Solution I-l 


Solutions II-(1,2) 




a = 1 


a = 2 


q;= 1 


a = 2 


M 


0.963t:«li 


l-630l;J°° 


0.965l:°lt 


1.00^_ Q4Q 




0.879l;l°° 


0.702* 


0.9011:1°° 


0.435* 


(a) 
94 


0.598t;J|° 


'-'■'-'-'--'--.060 


0.689t;Ji° 


0.687l:°o|° 




Pole position 


II sheet 


o.987i-;^f;-^ 

-/(o.()i5i;;;;;5) 


0.989l-"°5 


III sheet 


0.964l;«}^ 
-i(0.070l;°l°) 


1.5581;°^^ 
-i(0.14ll;°lg) 


0.9651:°!^ 
-i(0.073l;°l°) 


i.57ii:g^^ 
-i(o.i5il:°lg) 




ao-resonances with K-matrix background term 




Solution I-l 


Solutions II-(1,2) 




a = 1 


a = 2 


Q = 1 


a = 2 


M 


0.9441;°!^ 


1 624+"° 
-■-•o^^-.oso 


0.939l;°;«B 


l.D'±U_ Q4Q 




0.879i;l°° 


0.702* 


0.9011:1°° 


0.435* 




0.6511:^0 




0.653l:J«° 


0.65ll:°«° 




/ii = 0.52911°° 


so = 1.0+1° 


/n= 0.73111°° 


so = 1.9+1° 




Pole position 


II sheet 


0.990l:»°'^ 
-i(0.039l;°o«5) 


0.9931:°°^ 
-i(0.042l:°l°) 


III sheet 


0.965l:»l^ 
-i(0.063l;°l°) 


1.559t:°i 

-i(o.i45i:°i«5) 


0.965l:°l«5 
-i(0.068l:°l°) 


1.5751:°!;: 

-i(0.153l:°l^) 



23 



TABLE VII. 

Masses and coupling constants (in GeV) for /2 and 02 resonances. 





/2-resonances 




Solution I-l 


Solutions II-(1,2) 




a = l 


a = 2 


a = 1 


a = 2 


M 


1 236+°i° 

l.ZOO_ QJ^Q 


1.530l;«° 


1 oQQ+010 


i.529i:?,;g 




1.342l:i°° 


1.3421:1°° 


1.0381:}°° 


1.0381:}°° 




-(8.41^0) 


86.6l^J 


-(8.8ll;°) 


86.21^:^ 


9lJ 


0.318l:°^3S 




0.3181:;™ 


472+ ■°2f' 




a = TTTT 


a = rjri 


= TTTT 


a = rirj 


/la 


0.742l:°«,° 


-(i-oil;g^°o) 


0.2871:°™ 


-0.1431:}°° 


Ta 


1.997l;l^° 


1.077l;°g° 


2.4741:}^° 


1.295l:}fo 




/i3 = 0.684 ± 0.100 


/ba =0, & = 2, 3 So 


/i3 = 0.578 ± 0.100 

= 5.0 




Pole 
position 


1.262l:»g 
-i(0.092l;°°^) 


1.5141;'^;° 


1.2611:;;°;: 

-i(0.089l:°°^) 


-1 rr,r, + .005 

-z(0.076l:°°?) 




a2-resonances 




Solution I-l 


Solutions II- (1,2) 




a = 1 


a = 2 


a = 1 


a = 2 


M 


i.3i6l:°S? 


1.645l:°«,° 


1 019+. 005 
.'-•'J ."-^-.006 


1.6951:°^° 




i.o8ol:S! 


0.270l;l™; 


1.3001:12;; 


0.3251:};;;; 




0.38ll;°«,° 
ri = 1.8451:1^° 


0.5971;°^,° 


0.4261:°^° 
ri = 2.4061:}^° 


0.617l:°i° 


Pole 
position 


1.309l:°°«, 
-i(0.058i;°°^) 


1-6151:°^° 
-i(0.12li:°l^) 


1.3081:°°^ 
-i(0.059l:°°^) 


1.667l:g^g 
-i(0.123l:°}«) 
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